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^ ABSTRACT 

m 

<^ \ Cosmological implications of clusters of galaxies are discussed with particular attention to their impor- 
2^ ' tance in probing the cosmological parameters. More specifically we compute the number counts of clusters 
of galaxies, Log N - Log S relation, in X-ray and submm bands on the basis of the Press-Schechter theory. 
As an important step toward breaking the degeneracy among the viable cosmological models, we observed 
the most luminous X-ray cluster RXJ1347-1145 in three bands (21 and 43 GHz in the Nobeyama Radio 
Observatory, Japan, and 350 GHz in the J. C. Maxwell telescope at Mauna Kea, Hawaii). We report on 
^ \ the preliminary results which are in good agreement with the profile of the Sunayev - Zel'dovich effect 
f^ ' predicted on the basis of the previous X-ray observation of the cluster. 



ON 

Oh; 

6 



^ : 1 INTRODUCTION 

_ _ I 

There are several reasons why clusters of galaxies are regarded as useful tools in cosmology; (i) since the 
dynamical time-scale of clusters is comparable to the age of the universe, they should retain the cosmo- 
logical initial conditions fairly faithfully, (ii) clusters can be observed in various wavelengths including 
optical. X-ray, radio, mm and submm bands, and recent and future big projects (e.g., SDSS, AXAF, 
PLANCK) aim to make extensive surveys and detailed imaging/spectroscopic observations of clusters, 
(iii) to the first order approximation, clusters can be regarded as a simple system of dark matter, gas 
and galaxies, and thus theoretically well-defined and relatively well-understood, at least compared with 
galaxies themselves, (iv) clusters of galaxies can be observed up to high redshifts and thus provide a 
probe of the distant universe. 

X-ray observations are particularly suited for the study of clusters, because the X-ray emissivity is 
proportional to n^ and less sensitive to the projection contamination which has been known to be a 
serious problem in their identifications with the optical data. In addition, the Sunyaev - Zel'dovich 
(SZ) effect (Sunyaev &: Zel'dovich 1972), observed in radio, millimeter and submillimeter bands, is now 
opening a new window to investigating cluster properties, especially at high redshifts. In this paper, we 
aim to show the significance of current and future observations in these bands in cosmology. 



2 LOG N - LOG S OF X-RAY CLUSTERS 

We compute the number of clusters observed per unit solid angle with X-ray flux greater than S by 

Ni>S)= rdzd\{z)c f r dS {l + zfnMiM,z)^§^^, (1) 

where c is the speed of light, t is the cosmic time, dA is the angular diameter distance, Tgas and L^and are 
respectively the gas temperature and the band-limited absolute luminosity of clusters, and umIM, z)dM 
is the comoving number density of virialized clusters of mass M ^ M + dM at redshift z. 

Given the observed flux S in an X-ray energy band [Ea^E^], the source luminosity L^and at z in the 
corresponding band [Ea{l + z),Ei,{l + z)] is written as 

Lband[Ea{l + z),E,{l + z)] = A7rdl{z)S[Ea, E^], (2) 

where d^ = {l+z)'^dA is the luminosity distance. We adopt the observed L^oi—Tgas relation parameterized 
by 
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We take L44 = 2.9, a = 3.4 and (^ = as a fiducial set of parameters on the basis of recent observational 
indications (David et al. 1993; Mushotzky &; Scharf 1997). Then we translate Liyoi{Tgas) into the band- 
limited luminosity Lhand[Tgas, Ei, E2] by properly taking account of metal line emissions in addition to 
the thermal bremsstrahlung (we fix the abundance of intracluster gas as 0.3 times the solar value). 

Assuming that the intracluster gas is isothermal, its temperature Tgas is related to the total mass M by 
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where ks is the Boltzmann constant, G is the gravitational constant, nip is the proton mass, /i is the mean 
molecular weight (we adopt fi = 0.59), and 7 is a fudge factor of order unity which may be calibrated 
from hydrodynamical simulations or observations. We adopt 7 = 1.2 as canonical based on the recent 
simulation results (e.g. Bryan & Norman 1998). The virial radius rvh{M,Zf) of a cluster of mass M 
virialized at Zf is computed from A^ir, the ratio of the mean cluster density to the mean density of the 
universe at that epoch. We evaluate this quantity using the formulae for the spherical collapse model 
presented in Kitayama & Suto (1996b) and assuming for simplicity that Zf is equal to the epoch z at 
which the cluster is observed. Finally, we compute the mass function umIM, z)dM in equation (jl]) using 
the Press-Schechter theory (Press & Schechter 1974) assuming z = Zf as above. The effect oi Zf ^ z is 
discussed by Kitayama & Suto (1997) in this context, and the more general consideration oi zj ^ z is 
given in Lacey & Cole (1993), Sasaki (1994), and Kitayama & Suto (1996a, b). 



3 BREAKING THE DEGENERACY 

WITH THE SUNYAEV - ZEL'DOVICH EFFECT 

As has been realized for a while, a cosmological observation is often accounted for by a set of quite 
different models. In the present case, we found that the following models equally fit the Log N - Log S 
of ROSAT X-ray clusters (quoted numbers are la statistical and systematic errors, respectively): 

CT8 = (0.54 ± 0.02 ± 0.08) x j^-O-SS-o.soQo+o.ssQg ^^^ 



for Ao = 1 - Oo Cold Dark Matter (CDM) models, and 



as = (0.54 ± 0.02 ± 0.08) x j^-o.28-o.9iCo+o.680g 



(6) 



for Ao = CDM models. In order to exhibit this degeneracy among viable cosmological models, we 
select several examples listed in Table 1 which reproduce the ROSAT Log N - Log S in almost an 
indistinguishable manner (Fig.||a). 

Table 1: CDM models from the ROSAT X-ray Log N - Log S. 
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LOS 0.3 0.7 0.7 1.04 3.4 1.2 

O045 0.45 0.7 0.83 3.4 1.2 

El 1.0 0.5 0.56 3.4 1.2 

L037 0.3 0.7 0.7 0.90 3.4 1.5 

LOla 0.1 0.9 0.7 1.47 2.7 1.2 
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Figure 1: LogA^-logS relations of galaxy clusters for CDM models in (a) the soft X-ray (0.5-2 keV) band, 
(b) the hard X-ray (2-10 keV) band, and (c) the submm (0.85 mm) band. Lines represent the models 
listed in Table 1; L03 (solid), O045 (dotted). El (short dashed), LO37 (long dashed), and LOla (dot- 
dashed). Also shown in panel (a) are the la error bars from the ROSAT Deep Cluster Survey (Rosati 
et al. 1995, 1997), and the error box from the ROSAT Brightest Cluster Sample (Ebeling et al. 1997, 
1998). In panel (b), the number of clusters observed by HEAO 1 (Piccinotti et al. 1982) is indicated by 
the la error bar, and the counts of all X-ray sources detected by ASCA are plotted for reference as a 
histogram (Ueda et al. 1998) and downward arrows (Cagnoni et al. 1998). 



While this degeneracy can be broken by several methods, one possibility is to combine with observations 
of the Sunyaev & Zel'dovich (SZ) effect, the distortion of the cosmic microwave background (CMB) 
spectrum due to the inverse-Compton scattering in high temperature intracluster electrons (Sunyaev & 



Zel'dovich 1972, SZ). The specific intensity of tlie SZ signal at the rest frame of the cluster is given by 
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where /ip is the Planck constant, and y is the Compton y-parameter: 
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In the above expressions, T2.73 is the CMB temperature Tjq divided by 2.73K, x is the photon energy in 
units of kT^Q, and 
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By integrating equation (^) over a cluster, one obtains the excess flux from the cluster relative to the 
CMB: 
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Figure 2: The spectral dependence of the intensity change AI^/{iQy) induced by the SZ effect. The 
non-relativistic expression g(x) is shown as triangles, whereas the relativistic results are plotted as lines 
for T = 5 keV (sohd), 10 keV (dotted), and 15 keV (dashed). 



The spectrum of the SZ effect is completely specified by the function g{x) which is independent of cluster 
properties. Most importantly, g(x) is negative at x < 3.83 (z^ < 217GHz, A > 0.14cm), while positive at 
X > 3.83 (Fig.0). Clusters of galaxies therefore serve as negative (absorbing) sources at mm and cm bands. 



and as positive (emitting) sources at subnini band. At T ^ 10 keV, relativistic effects of electrons become 
important and the above expressions need to be modified. We thus apply the relativistic correction to 
the intensity change derived by Rephaeli (1995; see also Rephaeli &: Yankovitch 1997): At the observed 
wavelength of 0.85 mm, the above correction leads to 4%, 11% and 16% reduction in the excess flux, i.e., 
g{x), at T = 3 keV, 8 keV and 12 keV, respectively. While clusters with T ;^ 10 keV are rare and the 
above correction does not make signiflcant difference, we take it into account for completeness. 
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Figure 3: Contour maps on the rJo-cg plane in (a) spatially flat (Aq = 1 — f^o) CDM models, (b) open 
(Ao = 0) CDM models, (c) spatially flat CDM-like models with the flxed shape parameter (F = 0.25), 
and (d) open CDM-like models with F = 0.25. In all cases, h = 0.7 and the flducial set of intracluster 
gas parameters is assumed. Shaded regions represent the la significance contours from the soft X-ray 
(0.5-2 keV) logA^-logS". Dotted and solid lines indicate the contours of the number of clusters greater 
than S per steradian (10^, 10^, 10^ from bottom to top) with S = 10^^^ ergcm^^s^^ in the hard X-ray 
(2-10 keV) band and with Si, = 50 mJy in the submm (0.85 mm) band, respectively. Thick dashed lines 
represent the COBE 4 year results by Bunn & White (1997). 



Panels (b) and (c) in Figure |T| show the number counts of clusters of galaxies in the hard X-ray (2-10 
keV) and the submm (0.85 mm, due to the SZ effect) bands for models which reproduce the observed Log 
A^ - Log S in the soft X-ray band (Table 1). Clearly the predictions in submm bands are significantly 
different among the models. To see this clearly, we plot in Figure y contour maps of the cluster logA-logS" 
in different bands on the f^o ~ ^s plane. For the 0.5-2 keV band, the la significance level is plotted. 



while for the 2-10 keV and submm bands, the contours of the number of clusters per steradian (10^, 10^, 
10*^ from bottom to top) at S'(2-10 keV) = 10~^^ ergcm~^s~^ and £',^(0.85 mm) = 50 mJy are plotted, 
respectively. The a^ values derived from the COBE 4 year data (Bunn & White 1997) are also given for 
reference. Panels (a) and (b) show that the contours for the 2-10 keV band counts run almost parallel 
to the x^ contour of the 0.5-2 keV band counts. In this sense, future logA^-logS* data in the hard X-ray 
band will provide an independent consistency check of the ROSAT soft X-ray data. The shape of the 
submm logA^-log5 contours, on the other hand, is quite different, especially at high a^, and thus should 
place complementary constraints on f^o and a^. If the fixed shape parameter (F = 0.25) is adopted, 
the corresponding logA^-logS* contours are similar to the conventional CDM (F ~ Vtoh) case, except at 
r^o ~ 0.2 and fig ^ 0.8. By contrast, the COBE normalized ag is very sensitive to any changes in the 
spectral shape. 



4 SEARCH FOR THE SUBMM SZ EFFECT OF RXJ 1347-1195 

RXJ 1347.5-1145 {z = 0.45) is the most luminous X-ray cluster of galaxies known so far (Schindler et 
al. 1997). Its exceptionally high column density of electrons at the center implies that this is one of 
the most promising clusters for the detection of the first unambiguous SZ increment signal. This is 
clearly illustrated in Table 2 where we compute the expected central SZ intensity A/jy(0) at 0.85mm for 
several clusters with measured SZ temperature decrement in the Rayleigh - Jeans regime (except for our 
current target, RXJ 1347-1195). Wherever available, lo" errors in the parameters are shown which are 
also adopted to estimate the error in A/jy(0). 



Table 2: Expected central SZ intensity A/jy(0) at 0.85mm for several clusters 
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T [keV] 
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6c [arcsec] 


n0[10-24fcm-'^l 


AJy(0)[/i5,) mjy/boam] 


RXJ 1347-1145 


0.45 




0.56 ±0.04 


8.4 ±1.8 


9.4 


7.1 ±1.7 


CL0016+16 


0.555 


8 4+-^-^ 


0.81 


51 


0.65 


1 o+U-3 


A2163 


0.201 


14.6 ±0.5 


0.62 ±0.01 


72 ±3 


0.77 ±0.03 


4.5 ±0.3 


A2218 


0.171 


6.7 ±0.45 




60 ±12 


0.63 


1.2 ±0.3 


A2142 


0.0899 


8.68 ±0.12 


1.0 ±0.3 


221 ± 84 


0.493 ± 0.023 


1.6 ±0.8 


A478 


0.0881 


6.56 ±0.09 


0.67 ±0.03 


116 ±18 


0.675 ±0.124 


1.3 ±0.3 


A2256 


0.0581 


7.51 ±0.11 


0.80 ± 0.02 


320 ± 12 


0.251 ±0.013 


0.88 ± 0.06 


Coma 


0.0235 


9.10 ±0.40 


0.75 ± 0.03 


630 ± 36 


0.289 ± 0.004 


1.1±0.1 



On the basis of this analysis, we performed the search for the SZ effect towards the cluster RXJ1347- 
1145, at 21 and 43 GHz with NRO (Nobeyama Radio Observatory) 45m telescope between March 3 and 
15, 1998, and at 350 GHz with SCUBA installed on JCMT (J.C.Maxweh Telescope) 15m telescope on 
May 30 and 31, 1998. Although there have been a few claims of the detection of the SZ temperature 
increment for some clusters (e.g. Andreani et al. 1996; Lamarre et al. 1998), they measured the total flux 
alone which could be significantly contaminated by the dust in our Galaxy and also by submm sources 
in the field. Therefore it is essential to map the profile of a cluster so that one can test it against the 
predicted SZ signal profile. This is the reason why we attempted the observation with SCUBA which 
has a reasonable angular resolution, 15" in FWHM, with a FOV of 160". 

As a matter of fact, it turned out that a radio point source exists near the center of our target cluster. 
The flux of the central point source had been determined to be 47.6 it 0.45 mJy at 1.4GHz with VLA 
(Condon et al. 1988), and 11.0 it 0.7 mJy at 28.5GHz in the OVRO serendipitous survey of SZ effect 
(J.Carlstrom, private communication). We also observed the central source at 100 GHz with Nobeyama 
Millimeter Array (NMA) between May 19th to 21st, 1998, and found the flux to be 5.0 ± 1.5 mJy. The 



three data at z^ < 100 GHz are very well approximated by a single power-law: 

Fp{ij) = (56.1 ± 1.0)(i//lGHz)-°-49±002mJy, 



(12) 



where the quoted errors are la. We subtracted the contribution of the point source according to equation 
(12) from our data at 21 and 43 GHz. Although the point source might be a variable, the degree of 
variability is known to be fairly small for radio sources with steep spectrum (spectral index ;$ — 0.5; e.g. 
Eckart, Hummel & Witzel 1989). 
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Figure 4: Radial intensity profile towards RXJ1347 at (a) 21 and (b) 43 GHz observed at NRO. The 
estimated contribution of the central radio source (12.7mJy at 21 GHz, and 8.9 mJy at 43 GHz) is 
subtracted. Filled circles indicate our data with la error-bars. Solid lines indicate the prediction of the 
SZ signal (without the point source contribution) using the best-fit parameters in the X-ray observation. 



Figure H shows the profiles of the scans towards the cluster at 21 and 43 GHz after subtracting the point 
source. In both cases, a extended negative intensity profile is clearly visible. Especially at 21 GHz where 
the S/N is significantly higher than at 43GHz, the profile is quite consistent with the SZ profile expected 
from the X-ray observation. This is the first detection of the SZ temperature decrement towards the 
cluster RXJ 1347-1145. 

The observation of the cluster at 350 GHz with SCUBA was performed on May 30 and 31, 1998. The 
raw data were processed with REMSKY in SURF package to remove spatially correlated sky-noise. Since 
REMSKY systematically changes the base-level of the map up to the sky-noise level, the resulting zero- 
level is uncertain depending on the sky condition. We estimated la error of our base-level or DC offset, 
IdC) to be as large as 2.9 mJy/beam; the sky condition during our observation was not so good (the zenith 
optical depth at 350 GHz ranged around rsso = 0.46 — 0.60). Also we extracted 7 spurious contamination 
sources above a threshold of 3a using SExtractor package. The radially averaged profile of the cluster at 
350 GHz after subtracting the sources and the uncertain DC level is plotted in Figure |5|. Also plotted in 
the figure are the SZ intensity profiles at 350 GHz predicted from the latest X-ray observation with the 
point source of -Fp = and 3.2 mJy, as well as the point spread function of the 3.2 mJy source. Since 
normal radio sources often show a steep decline of flux at higher (submillimeter) , we regard the value, 3.2 



mJy, extrapolated from equation (12) as the upper limit of the flux. Again it is clear that the extended 
feature of the observed intensity profile cannot be described by the point spread function alone. 
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Figure 5: Radial intensity profile towards RXJ1347 at 350 GHz observed at JCMT/SCUBA. Filled circles 
indicate our data with Icr error-bars. Dotted curve shows the PSF of 3.2 mJy source, and dashed (solid) 
curves plot the SZ profiles with (without) the possible point source contribution, respectively, using the 
best-fit parameters in the X-ray observation. 
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Figure 6: Confidence contours on the y-parameter and the core radius 9c from combined data analysis 
at 21 and 350 GHz. The cross indicates the parameters determined from X-ray observations by ROSAT 
and ASCA satellites. 



To express the result in a statistical and quantitative manner, we carried out a y^ fit to the SZ profiles at 
21 and 350 GHz simultaneously on the basis of the /5-model with taking a central y-parameter, y^^ and the 
angular core radius, Q^^ as free parameters. The results are summarized in Figure ^. The agreement with 
the X-ray observation is remarkable and reassuring, and we conclude that the present data provide the 
strongest and most convincing case for the detection of the submm SZ signal from the cluster. Detailed 
analysis taking account of various astrophysical effects and the cosmological implications will be presented 
elsewhere (Komatsu et al. 1999a, b). 



5 DISCUSSION AND CONCLUSIONS 

It is widely recognized now that many cosmological models are more or less successful in reproducing 
the structure at redshift z ^ by construction. This is because the models have still several degrees of 
freedom or cosmological parameters which can be appropriately adjusted to the observations at z ~ {Qq, 
(Tg, h, Ao, b{r, z)). There are several ways to break such cosmological degeneracy. For instance, combining 
the COBE 4 year data with the cluster Log A^ - Log S observation, we showed that low-density CDM 
models with (f^o, Aq, /i, erg) = (0.3,0.7,0.7,1) and (0.45,0,0.7,0.8) are viable models, while the Einstein 
- de Sitter model is strongly disfavored even if equation (|5|) can be satisfied by adjusting cg accordingly. 

Along with such consideration, we proposed a cluster survey based on the SZ effect. As a first step toward 
the goal, we conducted the search for the SZ effect towards the cluster RXJ1347-1145 in multi-bands. 
We detected the extended signal which are in good agreement with the predicted SZ profiles, especially 
at 21 and 350 GHz. In particular, our results provide the strongest and most convincing case for the 
detection of the SZ temperature increment in the submm band in which the profile of a cluster of galaxies 
is resolved for the first time. 
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